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Abstract 

This paper is devoted to the extension to the full 3x3 Euler sys- 
tem of the basic analytical properties of the equations governing a fluid 
flowing in a duct with varying section. First, we consider the Cauchy 
problem for a pipeline consisting of 2 ducts joined at a junction. Then, 
this result is extended to more complex pipes. A key assumption in 
these theorems is the boundedness of the total variation of the pipe's 
section. We provide explicit examples to show that this bound is nec- 
essary. 

Keywords: Conservation Laws at Junctions, Coupling Conditions at 
Junctions. 
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1 Introduction 

We consider Euler equations for the evolution of a fluid flowing in a pipe 
with varying section a = a(x), see [T71 Section 8.1] or [P2| [T5]: 

d t (ap) + d x (aq) = 

d t {aq) + d x [aP (p, q, E)] = p (p, e) d x a (1.1) 
d t (aE)+d x [aF(p,q,E)] =0 

where, as usual, p is the fluid density, q is the linear momentum density and 
E is the total energy density. Moreover 

1 2 2 

E(p,q,E) = -^-+pe, P(p, q, E) = q ~+p, F(p,q,E) = -(E+p) , (1.2) 

2 p p p 

with e being the internal energy density, P the flow of the linear momentum 
density and F the flow of the energy density. The above equations express 
the conservation laws for the mass, momentum, and total energy of the fluid 
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through the pipe. Below, we will often refer to the standard case of the ideal 
gas, characterized by the relations 

p = (7 — l)pe, S = lne — (7 — 1) lap , (1-3) 

for a suitable 7 > 1. Note however, that this particular equation of state 
is necessary only in case (p) of Proposition 13.11 and has been used in the 
examples in Section HI In the rest of this work, the usual hypothesis [161 
formula (18.8)], that is p > 0, d T p(r,S) < and 8^ t p(t,S) > 0, are suffi- 
cient. 

The case of a sharp discontinuous change in the pipe's section due to 
a junction sited at, say, x = 0, corresponds to a(x) = a~ for x < and 
a(x) = a + for x > 0. Then, the motion of the fluid can be described by 

dtp + d x q = 
< d t q + d x P(p,q,E)=0 (1.4) 
d t E + d x F(p,q,E)=0, 

for x 7^ 0, together with a coupling condition at the junction of the form: 

* (a~, (p,q,E)(t,0-);a + ,(p, q,E) (t,0+)) =0. (1.5) 

Above, we require the existence of the traces at x = of (p, q, E). Various 
choices of the function \& are present in the literature, see for instance [U 
[U[9] in the case of the p-system and [10] for the full 3x3 system (jl.4p . Here, 
we consider the case of a general coupling condition which comprises all the 
cases found in the literature. Within this setting, we prove the well posedness 
of the Cauchy problem for f|1.4j> — f)1.5j> . Once this result is obtained, the 
extension to pipes with several junctions and to pipes with a W 1,1 section 
is achieved by the standard methods considered in the literature. For the 
analytical techniques to cope with networks having more complex geometry, 
we refer to [llj . 

The above statements are global in time and local in the space of the 
thermodynamic variables (p,q,E). Indeed, for any fixed (subsonic) state 
(p,q,E), there exists a bound on the total variation TV (a) of the pipe's 
section, such that all sections below this bound give rise to Cauchy problems 
for (|1.4p - (|1.5p that are well posed in L 1 . We show the necessity of this 
bound in the conditions found in the current literature. Indeed, we provide 
explicit examples showing that a wave can be arbitrarily amplified through 
consecutive interactions with the pipe walls, see Figure [TJ 

The paper is organized as follows. The next section is divided into three 
parts, the former one deals with a single junction and two pipes, then we 
consider n junctions and n + 1 pipes, the latter part presents the case of a 
W 1 ' 1 section. Section [3] is devoted to different specific choices of coupling 
conditions (|1.5|) . In Section [U an explicit example shows the necessity of the 
bound on the total variation of the pipe's section. All proofs are gathered 
in Section [5j 
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2 Basic Well Posedness Results 



Throughout, we let u = (p, q, E). We denote by M + the real halfhne [0, +00 [, 
while M + = ]0, +00 [. Following various results in the literature, such as [TJ 
[21 [H [ini 02], we limit the analysis in this paper to the subsonic region 
given by Xi(u) < < A3 (it) and \2(u) 7^ 0, where Aj is the i—th eigenvalue 
of (|1.4p . see ([5.ip . Without any loss of generality, we further restrict to 

A = |« G R+ x M+ x R+: Ai(it) < < A 2 (u)} . (2.6) 

Note that we fix a priori the sign of the fluid speed v, since A 2 (u) = q/p = 
v>0. 



2.1 A Junction and two Pipes 

We now give the definition of weak 'I'— solution to the Cauchy Problem 
for (|1.4p equipped with the condition (|l,5p . extending Definition 2.1] 
and Definition 2.2] to the 3x3 case (|1.4|) and comprising the particular 
case covered in [10\ Definition 2.4]. 

Definition 2.1 Let (M + x ^ ) 2 -> u Q G BV(M; A ) and two positive 
sections a~ , a + be given. A ^> -solution to ( fi.^| j with initial datum u Q is a 
map 

u e C°{R+;Ll oc (R+;A )) 
u(t) G BV(R;i ) for a.e. t£R+ K ' ' 

such that 



1. for x 7^ 0, u is a weak entropy solution to |j.^[ ); 

2. for a.e. 16K, u(0, x) = u Q (x); 

3. /or a.e. t £ M + , i/ie coupling condition hi. 5^ at the junction is met. 

Below, extending the 2x2 case of the p-system, see [H 19] , we consider 
some properties of the coupling condition (|1.5|) . which we rewrite here as 

^(a',u';a + ,u + ) = 0. (2.8) 
(*0) Regularity: * G C 1 x A,) 2 ; I 

('I'l) No-junction case: for all a > and all u~,u + G Aq, then 

^(a, a, u + ) = <^=^> u~ = u + . 

(\E r 2) Consistency: for all positive a~,a°,a + and all u~,u°,u + G Aq, 

^>(a~,u~;a°,u ) = T , _ _ + 

^(a°,n°;a + ,u + ) = v ' ' ' ' 
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Moreover, by an immediate extension of O Lemma 2.1], (\I/0) ensures that 
(|2.8p implicitly defines a map 

u+ = T(a-,a + ;u-) (2.9) 

in a neighborhood of any pair of subsonic states u~ , u + and sections a~ , a + 
that satisfy ^(a~, u~; a + , u + ) = 0. 

The technique in [Bj allows to prove the following well posedness result. 

Theorem 2.2 Assume that '9/ satisfies conditions (&())-(& 2). For every 
a > and u 6 ^4o site/t i/jai 

det ti(«) D u+ *-r 2 («) r 3 (ti)] ^0 (2.10) 

i/iere exist positive 5, L such that for all a - , a + with \a + — a| + |a _ — a| < <5 
i/iere exists a semigroup S: M + xP-»D w;zf/i i/te following properties: 

1. V^{ueu + L 1 ^; A ): TV(«) < o} . 

5. For a// u € D, Squ = u and for all t,s > 0, StS s u = S s +t u - 
3. For all u,u' £ V and for all t, t' > 0, 

\\Stu — Sfiu'W^ < L ■ (\\u — u'W^i + \t — t'^j 

4- If u € T> is piecewise constant, then for t small, Stu is the gluing of 
solutions to Riemann problems at the points of jump in u and at the 
junction at x = 0. 

5. For all u Q £ V, the orbit t — > Stu Q is a \P -solution to with initial 

datum u Q . 

The proof is postponed to Section [5j Above ri(u), with i = 1,2,3, are 
the right eigenvectors of Df(u), see (|5.ip . Moreover, by solution to the 
Riemann Problems at the points of jump we mean the usual Lax solution, 
see [31 Chapter 5], whereas for the definition of solution to the Riemann 
Problems at the junction we refer to [HJ Definition 2.1]. 

2.2 n Junctions and n + 1 Pipes 

The same procedure used in [91 Paragraph 2.2] allows now to construct the 
semigroup generated by (|1.4p in the case of a pipe with piecewise constant 
section 

n-l 

a = a X]-oo,xx] + a i Xfaw+ii + ° n X[x„,+oo[ 
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with n G N. In each segment ]xj, Xj+% [, the fluid is modeled by (|l,4p . At 
each junction Xj, we require condition (jl.5p . namely 

for all j = 1, . . . , n, where 
*(o i _ 1 ,uJ;a i ,«T) = n ± = li m±Ui ( x ). (2.11) 

We omit the formal definition of ^-solution to (|1.4p ~ (|1.5p in the present 
case, since it is an obvious iteration of Definition [2HJ The natural extension 
of Theorem 12.21 to the case of (|l,4p - ()2.1ip is the following result. 

Theorem 2.3 Assume that ^ satisfies conditions ( i SfO)-(ty2). For any a > 

and any u £ Aq, there exist positive M, A, 5, L, A4 such that for any pipe 's 
profile satisfying 

aePC(t;]ffl-A,o + A[) with TV (a) < M (2.12) 

there exists a piecewise constant stationary solution 

n-1 

U = UQXj-ooM + J2 ^X] XjtXj+1 [ + UnX]x n ,+oo[ 

i=i 

to ftl.4\ )- [2^TI\) satisfying 

Uj G Aq with \uj — u\ < 5 for j = 0, . . . n 
* [aj-\,Uj-i',a,j,Uj) = for j = 1,. . . ,n 
TV(ti) < M TV {a) (2.13) 

and a semigroup S a :M + x T> a — > T> a such that 

1. V a ^ {uefi + L 1 (R;Ao):TV(?i-'&) <6}. 

2. S3 is the identity and for all t,s>0, S?S% = S% +t . 

3. For all u, u' £ V a and for all t, t' > 0, 

||5fu-^u'|| Ll < L ■ (J|(u) -u'\\ Ll + \t-t'\\ . 

4- If u € T> a is piecewise constant, then for t small, Stu is the gluing of 
solutions to Riemann problems at the points of jump in u and at each 
junction Xj. 

5. For all u G V a , the orbit t —* Sfu is a weak ^-solution to |i.^| )-( fi07p . 

We omit the proof, since it is based on the natural extension to the present 
3x3 case of Theorem 2.4]. Remark that, as in that case, 5 and L depend 
on a only through a and TV (a). In particular, all the construction above is 
independent from the number of points of jump in a. 
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2.3 A Pipe with a W 1 1 Section 

In this paragraph, the pipe's section a is assumed to satisfy 

a G W 1 ' 1 (R; ]a - A, a + A[) for suitable A > 0, a > A 
TV (a) < M for a suitable M > 

a'(x) = for a.e. x G R \ [— X, X] for a suitable X > . 



(2.14) 



The same procedure used in [91 Theorem 2.8] allows to construct the semi- 
group generated by (jl.ip in the case of a pipe which satisfies (|2.14p . Indeed, 
thanks to Theorem 12.31 we approximate a with a piecewise constant func- 
tion a n . The corresponding problems to (|l,4p - (|2.1ip generate semigroups S n 
defined on domains characterized by uniform bounds on the total variation 
and that are uniformly Lipschitz in time. Here, uniform means also inde- 
pendent from the number of junctions. Therefore, we prove the pointwise 
convergence of the S n to a limit semigroup S, along the same lines in 
Theorem 2.8]. 



3 Coupling Conditions 

This section is devoted to different specific choices of ([2.8 



(S)-Solutions We consider first the coupling condition inherited from the 
smooth case. For smooth solutions and pipes' sections, system (jl.ip is equiv- 
alent to the 3x3 balance law 



dtp + d x q = — d x a 
a 



t 



d t q + d x P(p, q, E) =-—d x a 

d t E + d x F(p, q,E) = d x a. 

a 



The stationary solutions to (jl.ip are characterized as solutions to 



(3.1) 



d x (a(x)q)=0 (d x q 
d x (a(x) P(p,q,E)) = p(p,e)d x a or< 
d x (a(x)F{p,q,E)) =0 



4d r a 



d x P(p,q,E) 
d x F(p,q,E) 



fd x a 



(3.2) 



As in the 2x2 case of the p-system, the smoothness of the sections induces 
a unique choice for condition (|2,8p . see O (2.3) and (2.19)], which reads 



(S) * 



a + P{u + ) - a~P{u-) + / p(TZ a (x),£ a (x))a'(x)dx 

J-x 

a + F(u + ) - a~F{u~) 



(3.3) 
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where a = a(x) is a smooth monotone function satisfying a(—X) = a~ and 
a{X) = a + , for a suitable X > 0. 1Z a ,£ a are the p and e component in the 
solution to (|3.2p with initial datum u~ assigned at —X. Note that, by the 
particular form of (|3.3p . the function \& is independent both from the choice 
of X and from that of the map a, see 2. in Proposition 2.7]. 

(P)-Solutions The particular choice of the coupling condition in [10^ Sec- 
tion 3] can be recovered in the present setting. Indeed, conditions (M), (E) 
and (P) therein amount to the choice 



(P) ^(a~,u~,a + ,u + ) 



a q 

P(u-) 

a'F(u~) 



(3.4) 



where a + and a~ are the pipe's sections. Consider fluid flowing in a hori- 
zontal pipe with an elbow or kink, see [13] • Then, it is natural to assume 
the conservation of the total linear momentum along directions dependent 
upon the geometry of the elbow. As the angle of the elbow vanishes, one 
obtains the condition above, see [101 Proposition 2.6]. 

(L)-Solutions We can extend the construction in [H El [3] to the 3x3 
case (jl.4p . Indeed, the conservation of the mass and linear momentum in [3] 
with the conservation of the total energy for the third component lead to 
the choice 



(L) fy(a ,u ,a + ,u + ) 



a/q^ — a q 
a + P{u + ) - a"P(u 
a + F(u + ) - a~F(u- 



(3.5) 



where a + and a~ are the pipe's sections. The above is the most immedi- 
ate extension of the standard definition of Lax solution to the case of the 
Riemann problem at a junction. 



(p)-Solutions Following [H[2], motivated by the what happens at the hy- 
drostatic equilibrium, we consider a coupling condition with the conservation 
of the pressure p(p) in the second component of Thus 



(p) 9 (a ,u ,a + ,u + ) = 
where a + and a~ are the pipe's sections. 



cr q + — a q 
P(p + ,e + ) -p(p~,e~) 
a + F(u+) -a~F{u ) 



(3.6) 
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Proposition 3.1 For every a > and u G Aq, each of the coupling con- 
ditions * in < TOj) . ( fg^P , < TOj) . satisfies the requirements (*S> 0)-(*S> 2) 
and \2.10\) . In the case of \3. 6\) . we also require that the fluid is perfect, 
i.e. that U.3\) holds. 

The proof is postponed to Section^ Thus, Theorem 12 . 2 1 applies . yielding 
the well posedness of (|l,4p - (|1.5p with each of the particular choices of ^ 

in (HOD, (HID, dSSJ). 

4 Blow-Up of the Total Variation 

In the previous results a key role is played by the bound on the total variation 
TV(o) of the pipe's section. This requirement is intrinsic to problem (|1.4p 
(|1.5p and not due to the technique adopted above. Indeed, we show below 
that in each of the cases (|3.3p . (|3.4|) . (|3,5|) . (|3.6|) . it is possible to choose an 
initial datum and a section a E BV(1; [a~,a + ]) with a + — a~ arbitrarily 
small, such that the total variation of the corresponding solution to (jl.4p - 
(|1.5p becomes arbitrarily large. 

Consider the case in Figure [TJ A wave hits a junction where the 



a 




Aa 
















X 




Figure 1: A wave hits a junction where the pipe's section increases by 
Aa. From this interaction, the wave arises, which hits a second junction, 
where the pipe section decreases by Aa. 

pipe's section increases by, say, Aa > 0. The fastest wave arising from 
this interaction is , which hits the second junction where the section 
diminishes by Aa. 
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Solving the Riemann problem at the first interaction amounts to solve 
the system 




where u ^ Aq, see Figure [2] for the definitions of the waves' strengths af and 
cr^ . Above, T is the map defined in (|2.9p . which in turn depends from the 
specific condition f|2.8j) chosen. In the expansions below, we use the (p,q,e) 
variables, thus setting u = (p, q, e) throughout this section. Differently from 




Figure 2: Notation used in (|4"TTj) and (|P|) . 

the case of the 2x2 p-system in [9], here we need to consider the second 
order expansion in Aa = a + — a~ of the map T; that is 

T(a,a + Aa;u) = u + H(u)— + G{u) (— \ + o ( — ) (4.2) 

a \ a 1 V a I 



The explicit expressions of H and G in (|4.2p . for each of the coupling con- 
ditions (1231), (ELD, (1221), (GS3), are in Section O 

Inserting (|4.2p in the first order expansions in the wave's sizes of ()4.ip . 
with fi for i = 1, 2, 3 as in (|5.3p . we get a linear system in af, , ■ Now, 
introduce the fluid speed v = q/p and the adimensional parameter 

- \ 2 2 



7(7- l)e ' 



a sort of "Mach number". Obviously, i? G [0, 1] for u E Aq. We thus obtain 
an expression for of the form 

4 = f 1 + AW v + /2(,9) a * ■ (43) 



The explicit expressions of f\ and /2 in (|4.3p are in Section 15.21 

Remark that the present situation is different from that of the 2x2 
p-system considered in [9]. Indeed, for the p-system /2(#) = f2{'& + ) = 0, 
while here it is necessary to compute the second order term in (Aa)/a. 
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Concerning the second junction, similarly, we introduce the parameter 
= (t> + /c + ) 2 which corresponds to the state u + . Recall that u + is defined 

by u + = L3 \T (l%(u; crj); j , see Figure [2] and Section [531 for the 

explicit expressions of $ + . We thus obtain the estimate 



(7 



3 



l-/i(^ + )^ + / 2 (^ + ) (^) \t, (4-4) 



where i9 + = -!? + ^1?, a 3 , (Aa)/aj. Now, at the second order in (Aa)/a and 
at the first order in , (|4.3|) and (|4.4|) give 



( Wl( ^ + /2( , +) (^) 2 ) 

*(i + Aw£ + M*> 
= fi + x(*)(^)V 3 -- («) 

Indeed, computations show that /1 (1?) — /1 (^ + ) vanishes at the first order 
in (Aa)/a, as in the case of the p-system. The explicit expressions of x are 
in Section 15.21 

It is now sufficient to compute the sign of x- If h is positive, then 
repeating the interaction in Figure [T] a sufficient number of times leads to 
an arbitrarily high value of the refracted wave 03 and, hence, of the total 
variation of the solution u. 

Below, Section [5] is devoted to the computations of x m the different 
cases (|3.3p , (|3.4p , (|3.5p and (|3.6p . To reduce the formal complexities of the 
explicit computations below, we consider the standard case of an ideal gas 
characterized by (|1.3p with 7 = 5/3. 

The results of these computations are in Figure [3l They show that in 
all the conditions (|1.5p considered, there exists a state u S Aq such that 
> 0, showing the necessity of condition (|2.12p . However, in case (L), 
it turns out that x is negative on an non trivial interval of values of 1?. If 
u is chosen in this interval, the wave 03 in the construction above is not 
magnified by the consecutive interactions. The computations leading to the 
diagrams in Figure [3] are deferred to Section 15.21 
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1000 



10000 




theta in case (L) theta in case (p) 

Figure 3: Plots of x as a function of •&. Top, left, case (S); right, case (P); 
bottom, left, case (L); right, case (p). Note that in all four cases, x attains 
strictly positive values, showing the necessity of the requirement (|2.12p . 
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5 Technical Details 



We recall here basic properties of the Euler equations (jl.ip . (jl.4p . The 
characteristic speeds and the right eigenvectors have the expressions 



A, 



A 2 





-p 




P 




P 


n = 


pc-q 


T2 = 


q 


rs = 


q + pc 




qc — E — p 








E +p + qc 



(5.1) 



whose directions are chosen so that VAj • r j > for i = 1,2, 3. In the case 
of an ideal gas, the sound speed c = d p p + p~ 2 p d e p becomes 

c= V7(7-l)e- (5-2) 
The shock and rarefaction curves curves of the first and third family are: 



Si(u ,a) 



-a + Po 



V = V - 

e = e -\(p + p ) 
o- + Po 



for 



l i_ 

P Po 



P 

v = v 



(P ~ Po) 



1 1_ 

P Po 



for 



e = c ~ \ (p + Po) 



l i_ 

P Po 



Ri(u ,a) 



P = ~ 

V = V, 



V + Po 
•P 



[(pc)(p,S )]-'dp 



for 



Po 



R 3 (u ,a) 



S(p,e) = S(p ,e ) 
p = a + po 



v = v Q + 



[(pc)(p,S )]- L dp 



for 



S(p,e) = S(p ,e ) 
The 1,2,3-Lax curves have the expressions 



a 


< 





P 





Po 


V 


< 


V 


S 


> 


So 


a 


< 





P 


< 


Po 


V 


< 


Vo 


S 


< 


So 


a 


> 





P 





Po 


V 


> 




e 


< 


e 


a 


> 





P 


> 


Po 


V 


> 


Vo 


e 


> 


e 



Pot Qoi E ) 
L 2 (a;p ,q ,E ) 
L 3 (a;p ,q ,E ) 



Sx(a;p ,q ,E ), a < 
Ri(a;p 0l q ,E ), a > 

p = a + po 

v = v Q 

p(p, e) = p{p , e Q ) 

Sz{p;pi,qi,Ei), a < 
R 3 (a;p 0l q ,E ), a > 
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Their reversed counterparts are 



S 1 (u ,a) 



p = a + p 

V = V + J-(p- Po ) (i-i 



for 



e = e. 



o-Up + Po) (j-i) 



S3 (u Q ,a) 



-<? + Po 



v = v + 



(P - Po) 



1 1_ 

P Po 



for 



e = e D - \ (p + po) 



1 i_ 

P Po 



Ri (u ,<j) 



p = a + po 



V = V n 



[(pc){p,S )} dp for 



Po 



-R 3 ( M o,0") 



S(p,e) = S(p ,e ) 

p= -cr + po 
r-p 

v = v a + [(/?c)(p, S'o)] -1 dp 

-'Po 

5(/9,e) = S(p ,e ) 



for 



(7 




U 


P 


< 




u 


> 


l> 


S 


< 




(7 


< 





P 




Po 


?; 





Vo 


S 


> 


So 


a 


> 





P 


> 




V 


< 




e 


> 




a 


> 





P 


< 


Po 


V 


< 




e 


< 





and 



Li(a;p ,q ,E ) 

^2 (o-; p ,q ,E ) 

L 3 (cr; p ,q ,E ) 



S t (a; p ,q ,Eo), a < 

Poi 9o) 
P = -O + Po 

■y = v 

p{p,e) =p(p ,e ) 
Sz{<J\Po,q ,E ), a < 

R 3 {<y,Po,qo,E ), cr>o. 



In the (/9, g, e) space, for a perfect ideal gas, the tangent vectors to the 
Lax curves are: 



1 



- ^/7(7- 1)< 
-(7-l)f 





1 " 
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, f 2 = 


P 


, ^3 = 




e 






P. 





1 



J - ^7(7 " I> 
(7-l)f 



. (5.3) 



5.1 Proofs of Section [2] 

The following result will be of use in the proof of Proposition 
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Proposition 5.1 Let Oi \—* Li(uQ,ai) be the i-th Lax curve and Oi \—* 
L^~(ito,<Jj) be the reversed i-th Lax curve through uq, for i = 1,2,3. The 
following equalities hold: 



din 

do-i \ ai=0 



( 



1 

\\{u ) 

E D + Po q 



\ 



dL, 



' da 



\ Po Po ) 

I 



2 | (72=0 



dLs 



da. 



3 | o-3=0 



1 

A 3 (if ) 
E +Po , Qo 



( 1 

A2(n G ) 

Eq+P q Po Co 
V Po 

\ 



\ 



dePo J 



\ +—C , 
\ Po Po J 



+ 







1 

dLi 



dp | CTi =o dp o \ a . =0 dq a | CT . =0 dq o \ a . =0 dE o \ a . =0 dE o \ a . =0 
The proof is immediate and, hence, omitted. 

Proof of Theorem \2."A Following [7, Proposition 4.2], the 3x3 sys- 
tem (|1.4p defined for x E M. can be rewritten as the following 6x6 system 
defined for x E M + : 



d t U + d x T(U)=0 (t,x) E 
b (U(t, 0+)) = t E 



(5.4) 



the relations between U and ii = (p,q,E), between J- and the flow in (jl.4[) 
being 



C/(t,x) 



p(t,-x) 
q{t,-x) 
E(t, -x) 
p(t,x) 
q(t,x) 
E(t,x) 



and J"(C7) 



#2 

P{UuV 2 ,Us) 
F(U 1 ,U 2 ,U 3 ) 
U 5 

P(U 4 ,U 5 ,U 6 ) 
F(U A ,U 5 ,U 6 ) 



with x E M + and E,P,F defined in (|1.2p ; whereas the boundary condition 
in Q53J) is related to (fT3|) by 

b{U) = * (o~ (C/i, C/ 2 , 17; 3 ); a+, (f/ 4 , C/ 5 , tf 6 )) 
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for fixed sections a~ and a + . 

The thesis now follows from [6, Theorem 2.2]. Indeed, the assump- 
tions (7), (b) and (f) therein are here satisfied. More precisely, condi- 
tion (7) follows from the choice (|2.6p of the subsonic region Aq. Simple 
computations show that condition (b) reduces to 



det 



det 
det 



d<J\ | CT1= o u da 2 1 0-2=0 u 9u+ 1 0-3=0 9a 3 |o 3 =0 

D u - * • n (u) - D u + * • r 2 (u) - D u + * • rs («)] 
D u - * • n (u) D u + * ■ r 2 (u) £> u + * • r 3 (it)] , 



which is non zero for assumption if u E Aq and a > 0. Condition (f) needs 
more care. Indeed, system (|5.4p is noi hyperbolic, for it is obtained gluing 
two copies of the Euler equations (|1.4p . Nevertheless, the two systems are 
coupled only through the boundary condition, hence the whole wave front 
tracking procedure in the proof of jSJ Theorem 2.2] applies, see also [3 
Proposition 4.5]. □ 

Proof of Proposition lff.il It is immediate to check that each of the 
coupling conditions (|3.3p . (|3.4p . (|3,5p . (|3,6p satisfies the requirements (*0) 
and (*1). 

To prove that (*2) is satisfied, we use an ad hoc argument for condi- 
tion (S). In all the other cases, note that the function 'I' admits the rep- 
resentation ^(a _ , u~; a + , u + ) = ip(a~ , u~) — ip(a + , u + ). Therefore, (\E r 2) 
trivially holds. 

We prove below (|2.10p in each case separately. Note however that for 
any of the considered choices of 

D u+ ^(a~,u~,a + ,u + )\ u= ^ a=a = -D u -^(a~, u~, a+, u + )\ u= ^ a=n (5.5) 

so that (|2.10p reduces to 

det [D u - * • n («) D u+ * • r 2 (u) D u + V ■ r 3 (u)] 
= — det D u +^ ■ det [ri(u) r 2 (u) ^(n)] . 

Thus, it is sufficient to prove that det D u +^(a~, u~, a + , u + )\ u=u ^ a=a / 0. 

(S)-solutions To prove that the coupling condition (|3.3p satisfies (^2), 
simply use the additivity of the integral and the uniqueness of the solution 
to the Cauchy problem for the ordinary differential equation 
Next, we have 

D u ( f P (K a (x),£ a (x)) a'{x) dx] = C 



'-X 



\u=u, a=a 
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since a'(x) = for all x, because a = a + = a. Thus, \P in (|3.3p satisfies 



D u +^f((l ,U , (2. , U ) \u=u, a=a 





a 3 det 



8eP | 






P > 


J? 


-1) 


a 2 




_ e +p 




I) 


P 



deP 



12 
p V p 



9 e p <? 2 
P I 2 " 





9eP 



| I 9e_P | } 



P V P 



= -a 3 Ai(-u.) A 2 (n) A 3 (u), 
which is non zero if u G Aq. 

(P)-solutions Concerning condition (|3.4p . we have 

L» u +^(a-,u",a + ,u + )| u=Sja=s 

a 



det 



£ f 2 — 
p V p 





9e£ 



of + 



9 e p / q 



P VP 2 " 



E+p 
P 



p p p 2 



;-E+P 



-a|(f + 1 



= -a 2 Ai(M)A 2 (n)A 3 (M), 
which is non zero if ft G Aq. 



(L)-solution For condition 
above case: 



the computations very similar to the 



D u +V(a ,u ,a + ,u + )\ u= ^ a= a = -a 3 X 1 (u)X 2 (u)X 3 (u), 
which is non zero if u G Aq. 

(p)-solution Finally, concerning condition (|3,6p . 

D u +^(a~ ,u' ,a + ,U + )\ u=Qya= a 

a 



det 



p / p p p 2 





Sep 
P 

p V p 



a 2 A 2 (n) ( C 2 + A|(n)^ 



which is non zero if u G Aq and if the fluid is perfect, i.e. ()1.3|) holds. □ 
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5.2 Computation of x m (14.51) 

The Case of Condition (S) Let ^ be defined in (|3.3p and set 

r x 

E(a',a + ,u) = / P {Tl a {x),S a {x)) a (x)dx 
J-x 

where the functions lZ a ,£ a have the same meaning as in (|3.3p . A pertur- 
bative method allows to compute the solution to f|3.2|) with a second order 
accuracy in (Aa)/a. Then, long elementary computations allow to get ex- 
plicitly the terms H and G in (|4.2p of the second order expansion of T: 



H(p,q,e) 



G(p,q,e) 



tf3_4tf 2 + 5i?-2 
#-3tf 2 + 3tf-l 9 

q 

2(-tf 3 + 2tf 2 -i?) 
3(tf 3 -3tf 2 + 3tf-l) 6 

4(i? 3 - 2 1? 2 ) 
3(tf 3 -3tf 2 + 3i?-l) P 

70 1? 4 - 257 1? 3 + 342 2 - 207 ■& + 36 



18 (tf 3 -3tf 2 + 3tf-l) 
Moreover, the coefficients /i, /2 in (|4.3p read 
-30+ (0-3) 



6V^(i?-l)-6(^-l) 
\A? (126 tf 4 - 505 1? 3 + 758 tf 2 - 489 ■& + 270) 



72 (y$ (tf 3 -3tf 2 + 3tf-l) -tf 3 + 3tf 2 -3tf + l 
42i? 4 - 183 1? 3 + 278tf 2 + 33^ + 54 



+- 



72 (V# (?? 3 -3i9 2 + 30-l) -0 3 + 30 2 -30 + l 



Next, x is given by 

(l26 19 4 - 506 19 3 + 773 19 2 - 480 1? + 279^ + 42 tf 4 - 174 1? 3 + 311 tf 2 + 96 ■& + 45 



X = 



36 (Vd (i? 3 - 3 tf 2 + 3 1? - 1) - i9 s + 3 tf 2 - 3 + l) 



The Case of Condition (P) Let ^ be defined in (|3.4p . With reference 
to (j4.2|) . we show below explicitly the terms H and G in (|4.2p of the second 
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order expansion of T, 



H(p,q,e) 



G{p,q,e) 



8(-0 3 + 20 2 -0) 
3(0 3 -30 2 + 30-l) P 

2(50 4 - 70 3 -0 2 + 30) 
9(0 3 -30 2 + 30-l) 6 

64 (0 3 + 30 2 ) 
27 (0 3 -30 2 + 30-l) P 

565 4 - 1599 1? 3 + 927 2 - 405 



81 (0 3 - 30 2 + 30 -l) 
Moreover, the coefficients /i, /2 in (|4.3p read 

yfd (9 2 + 2 - 27) + 3 2 - 42 - 9 



AW 



18\/0 (0- 1) - 18(0-1) 
V# (154 5 + 931 1? 4 - 4416 3 + 6570 2 + 990 + 89l) 

324 (\/0 (0 3 -30 2 + 30-l) -0 3 + 30 2 — 3 ?? + l) 

86 5 - 311 4 - 752 3 + 7038 2 + 1026 + 81 
324 (^0 (0 3 -3tf 2 + 30-l) -0 3 + 30 2 -30 + l) 

Next, x is given by 

\/0 (407 5 + 1931 4 - 7858 3 + 14766 2 + 1179 + 1863) 

X = 



+- 



The Case of Condition (L) Let * be defined in (j3"3jl . Then, 



40 



324 fV0 (0 3 -30 2 + 30-l) -0 3 + 30 2 -30 + l 
-23 5 + 141 4 + 2002 3 + 15714 2 + 2565 + 81 
324 (a/0 (0 3 -30 2 + 30-l) -0 3 + 30 2 -30 + l 



H(p,q,e) 



and G(p, Q, e) 



3(0- 1) r 

35 2 - 9(40-l) 
9(0-1) 



The coefficients /i, /2 in (j4.3j) read 

AW = o 

yfd (63 2 - 106 + 27) +210 2 -780 + 9 



AO?) 



36 



1) -0 + 1 
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so that x is 



X 



V$ (63tf 2 - 106^ + 27) + 21 t? 2 -78tf + 9 
18 (V# (0 



The Case of Condition (p) Let ^ be defined in (|3.6p . With reference 

to flMD, 



H(p,q,e) 



G(p,q,e) 



2(4tf 3 + 12tf 2 + 9tf) 
~4(2tf 3 + 9tf 2 ) +27(20 + 1) 

2(40 3 + 12 2 + 90) 
4(20 3 + 9tf 2 ) + 27(2tf + 1) 6 

4(0 3 + 3i? 2 ) 
~4(20 3 + 90 2 ) +27(20 + 1) 

12 (0 3 + 2 2 ) 
4(20 3 + 9tf 2 ) + 27(2tf + 1) 6 



with /i and fa given by 

-20 2 + 4t?l + 30-9 



so that 

x -- 



2 (40 2 + 120 + 9) 
32 4 + 8\/0 (40 3 + 90 2 -90) +316 1? 3 + 558 2 + 2160 + 81 
6 (16 4 + 96 3 + 216 2 + 216 + 81) 

m^+VQyfti (0 3 +0 2 - 30) + 7OO0 3 + llO70 2 -540 + 81 
6 (160 4 + 960 3 + 2160 2 + 2160 + 8l) ' 
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sf& ((0-3+6) i9 2 + 18i9-9a ; 7^ + 2(3-20-3") tf 2 + 6 (2 o^ + 3) 



# A 

Aa 



9v^ (1? - 1) + 9 (1? - 1) a 
\A? (l4 (ll ovT -30^ tf 5 + (990 - 301 0-3- ) i? 4 + 3^25a ; 7 - 236) tf 3 + (ill ovT -18) i9 2 + 45 (^12 - 0-3") tf- 378 0-3"^ 



2 

) 

a 



108 (v^? (i? 3 - 3 19 2 + 3 1? - 1) + i? a - 3 19 2 + 3 1? - 
4 (217(7,7 - 105^ ■a 5 + 2 (495 - 1489ctj) # 4 + 4 (928 - 177) tf 3 - 2 (l077 o\7 -9) tf 2 + 24 (260-3" + 15) J? /Aa\ 2 

108 (^ (tf 3 -3tf 2 + 3i7-l) + tf 3 -3i7 2 + 3tf-l) U/ 

((1I0-3" -30) tf 3 + (-19^ - 108) i9 2 +9 (50-3" - 6) j? + 270-3-^ +2 (310-3" - 15) j? 3 - 36 (ovT +3) # 2 - 18 (5^ +3) 



$ A 
Aa 



27\/tf (1? - 1) + 27 (tf - 1) 



\f& (2 (2330-3- - 300) t? 6 + (12790-3- - 5310) 1? 5 + (5400 - 2543 0-3") tf 4 + 6 (677o-- + 1242) # 3 + 36 (l09 - 2970-3-) i9 2 + 729 (2 - 5o\r) ■& - 1215 a 3 



486 I 



(V& (tf 3 -3tf 2 + 3tf-l)+tf 3 -3i7 2 + 3tf-l) 
g 4 (601 0-3- - 150) j? 6 + 2 (3521 0-3- - 2655) i9 5 + 8 (225 - 3814 o^) i9 4 + 36 (655 0-3" + 207) # 3 + 108 (530-3- +36) tf 2 + 162 (9a^ +25) ■& f AaV 



Aa 



486 (\/tf (tf 3 -3i7 2 + 3tf-l)+tf 3 -3tf 2 + 3tf-i) 

(^(770-3- - 210) tf 3 + (-130-3- -36) i9 2 + 27 (0-3- +2) 27o- ; 7^ 2+ (2170-3" - 105 ) ^ 3 - 4 (147a,- + 9) tf 2 + 18 (50-3- +3) ■& 

54 %/# (1? - 1) + 54 (■& - 1) \ij 



2 



i?+ = 1? + 



-2 ((7,7+6) # 3 + \/# (lO 0-3" tf 2 + 27 0-3" ■& + 27 a.r) +3(0-3- -18) tf 2 - 9 (a~ + 6) 1? Afl 

3 (4 19 2 + 12 15> + 9) a~ 
-48 (0-3- +5) tf 5 + (1440-3- ^ 4 + 780 0-3- tf 3 + 22140-3- + 1944o-3-i9 + 1215o-3-) / Aa \ 2 

48 (2i9 4 + 12tf 3 + 27 1? 2 + 27 1? + 12) \~~a~ 
-4 (91 0-3- + 342) i9 4 - 54 (ll 0-3- +56) i9 3 -216 (20-3" + 15) i5» 2 - 81 (50-3- + 18) / A 



, 2 

a 



48(2tf 4 + 12# 5 + 27tf 2 + 27#+12) \ a 

Above are the values of i? + in the cases (S), (P), (L) and (p). 
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